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TECHNICAL REPORT 



The RICIS Concept 


The University of Houston-Clear Lake established the Research Institute for 
Computing and Information Systems (RICIS) in 1 986 to encourage the NASA 
Johnson Space Center (JSC) and local industry to actively support research 
in the computing and information sciences. As part of this endeavor, UHCL 
proposed a partnership with JSC to jointly define and manage an integrated 
program of research in advanced data processing technology needed for JSC’s 
main missions, including administrative, engineering and science responsi- 
bilities. JSC agreed and entered into a continuing cooperative agreement 
with UHCL beginning in May 1986, to jointJy plan and execute such research 
through RICIS. Additionally, under Cooperative Agreement NCC 9-16, 
computing and educational facilities are shared by the two institutions to 
conduct the research. 

The UHCL/RiCIS mission is to conduct, coordinate, and disseminate research 
and professional level education in computing and information systems to 
serve the needs of the government, industry, community and academia. 
RICIS combines resources of UHCL and its gateway affiliates to research and 
develop materials, prototypes and publications on topics of mutual interest 
to its sponsors and researchers. Within UHCL, the mission is being 
implemented through interdisciplinary involvement of faculty and students 
from each of the four schools: Business and Public Administration, Educa- 
tion, Human Sciences and Humanities, and Natural and Applied Sciences. 
RICIS also collaborates with industry in a companion program. This program 
is focused on serving the research and advanced development, needs of 
industry. 

Moreover, UHCL established relationships with other universities and re- 
search organizations, having common research interests, to provide addi- 
tional sources of expertise to conduct needed research. For example, UHCL 
has entered into a special partnership with Texas A&M University to help 
oversee RJCIS research and education programs, while other research 
organizations are involved via the “gateway* concept 

A major role of RICIS then is to find the best match of sponsors, researchers 
and research objectives to advance knowledge in the computing and informa- 
tion sciences. RICIS, working Jointly with its sponsors, advises on research 
needs, recommends principals for conducting the research, provides tech- 
nical and administrative support to coordinate the research and integrates 
technical results into the goaJs of UHCL, NASA/JSC and industry. 
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Summary 


This final report documents the first-year research efforts for developing 
on-off pulse control techniques for flexible space vehicles which are 
sometimes required to maneuevr as quickly as possible with a minimum of 
structural vibrations during and/or a maneuver, particularly, in the presence 
of structural mode uncertainty. 

The study objective was to explore the feasibility of computing open-loop, 
on-off pulse control logic for uncertain flexible spacecraft. The results 
indicate that the proposed robustification or desensitization approach does 
generate robust on-off pulse sequences for uncertain flexible spacecraft. 

On the contrary to a common notion, the results show that properly 
coordinated, on-off pulse sequences can achieve a fast maneuvering time with 
a minimum of structural vibrations during and/or after a maneuver, even in 
the face of spacecraft modeling uncertainty. The time-optimal responses have 
been desensitized at the expense of the increased maneuvering time. 
However, it is emphasized that simply prolonging the maneuver time does 
not help to reduce residual structural vibrations caused by modeling 
uncertainty; a proper coordination of pulse sequences is necessary. 
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Abstract 

Both feedforward and feedback control approaches for uncertain dynamical sys- 
tems (in particular, with uncertainty in structural mode frequency) are investigated. 
The control objective is to achieve a fast settling time (high performance) and ro- 
bustness (insensitivity) to plant uncertainty. Preshaping of an ideal, time-optimal 
control input using a tapped-delay filter is shown to provide a fast settling time with 
robust performance. A robust, non-minimum-phase feedback controller is synthe- 
sized with particular emphasis on its proper implementation for a non-zero set-point 
control problem. It is shown that a properly designed, feedback controller performs 
well, as compared with a time-optimal open-loop controller with special preshaping 
for performance robustness. 
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1. Introduction 


Flexible structures, including robot manipulators and optical pointing systems in 
space, are sometimes required to reorient or reposition as quickly as possible with a 
minimum of residual structural vibrations. The control task for such systems becomes 
more difficult if they have many flexible modes within the rigid-body control band- 
width. The rapid maneuvering control problem of flexible dynamical systems has 
been investigated by many researchers, and various feedforward /feed back approaches 
for minimizing residual structural vibrations have been developed (e.g, see [1-9]). 

The basic idea behind the various feedforward approaches is to find an input forc- 
ing function (e.g, such as a versine function) which begins and ends with zero slope 
so that structural modes are less likely to be excited. Such an input function, how- 
ever, does not fully utilize the available maximum maneuvering force, and results 
in a slower response time and also residual structural vibrations. Most feedforward 
approaches (in particular, an open-loop optimal approach) require accurate modeling 
of the system and thus are not robust to plant modeling uncertainty. Recently, Singer 
and Seering [7-8] have developed an alternative, robust approach of shaping a feedfor- 
ward input command by acausally filtering out the frequency components of flexible 
mode resonances. Some tradeoffs, however, must be made between performance (a 
fast settling time) and robustness (insensitivity) to plant parameter uncertainty even 
for this preshaped open-loop approach. 

This paper provides a comparison between a preshaped feedforward command 
generator [7-8] and a robustified feedback controller [10-13] with non-zero set-point 
command for a reference-input tracking problem. It is, however, emphasized that one 
of the primary motivations for the use of closed-loop rather than open-loop control 
systems in practice is to cope with unexpected disturbances, which an open-loop 
controller cannot. Such a disturbance rejection problem is not considered in this 
paper. A simple generic model of uncertain dynamical systems shown in Fig. 1 is 
used to illustrate the control concepts and methodologies. 

The remainder of this paper is organized as follows. In Section 2, time-optimal 
control inputs are determined for this example problem. In Section 3, such time- 
optimal control inputs are preshaped using a tapped-delay filter to provide a rapid 
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maneuver and robust suppression of residual structural vibrations. A robust Hoc 
compensator design is discussed in Section 4, with special emphasis on a proper 
implementation of a non-minimum-phase compensator for a non-zero set-point control 
problem. It will be shown that a properly designed, feedback controller performs 
well, as compared with a time-optimal open-loop controller with special preshaping 
for parameter robustness. 


2. Time-Optimal Control 

Consider a generic example of flexible dynamical systems as shown in Fig. 1. It 
is assumed that the two bodies both have nominal unit mass (m! = m 2 = 1) and 
are connected by a spring with nominal stiffness k = 1. A control input force u acts 
on body 1 and the position of body 2 is the output to be controlled (y = x 2 ). It is 
assumed that the control input is bounded as |u| < 1. 

In this section, time-optimal, open-loop control input u(<) is determined by min- 
imizing the performance index 

J=[ t, dt = t j . (1) 

w’here tj is the final time to be obtained. The resulting time-optimal control input 
will be preshaped in Section 3, using a tapped-delay filter, to improve peformance 
robustness with respect to plant parameter uncertainty. 

Rigid-Body Time-Optimal Control 

For a “rigidized” model of the nominal system shown in Fig. 1, the equation of 
motion is simply 

(mi+m 2 )y=u. (2) 

The rest-to-rest, time-optimal solution for y( 0) = 0 and y{tj) = 1 can be found as 

u(() = u s (t) - 2 u s (t - !+) + i i,(t - t } ) (3) 

where tj = 2y(mi + rn 2 )y(tj) = 2.828 sec and u,(<) represents a unit-step function. 
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If this time-optimal input force is exerted on the nominal system with a flexible 
mode (Fig. 1), a significant residual structural vibration will occur, as can be seen in 

Fig. 2. 

Flexible-Body Time-Optimal Control 

Consider a time-optimal control problem for the flexible-body model shown in 
Fig. 1. The equations of motion are 

m\i\ -f k(x l — x 2 ) = u (4a) 

m 2 x 2 - k(x 1 - x 2 ) = 0 (4b) 


can 


where x, and x 2 are the positions of body 1 and boc 
also be represented in state-space form as 

x(t) = Ax(t) + Bu(t) 
l i{t) = Cx(t) 


(5a) 

(5b) 


where 


x = 
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X| X2 X 3 X4 
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C = [ 0 1 0 0 ] . 

The time-optimal control input for this problem, which first appeared in [10], can 
be solved using the following approach discussed in [4]. For our two-mass-spring 
model, three switching times are required. The nominal system with m l = m 2 = k = 
1 is first transformed into modal equations by the coordinate transformation 

( 6 ) 

where <7, and q 2 are the modal coordinates, and the resulting modal equations are 

= u/2 ( 7a ) 




' 1 

1 


9 1 

x 2 [ 


1 

-1 


. 92 


4 


<72 + = u / 2 


(7b) 

where a? = y/2 rad/sec is the nominal flexible mode frequency. For given boundary 
conditions 

?i(0) = 0, <71 (0) = 0, 

92(0) = 0, <72(0) = 0, 

9i(</) = 1, 9i(f/) = 0, 

= 0, <72 (tf) = 0, 

the following two constraint equations can be obtained as 

1 — 2cosu}t\ + 2cosut 2 — 2cosut 3 + cosutj = 0 
4 - t) + 2 (t f - Uf - 2 (tj - t 2 f + 2 (t f - t 3 ) 2 = 0 

where t ly t 2 , and t 3 are the switching times and tj the final time to be solved. From 
the symmetric nature of this problem, we have 

t 2 = tj/2 

— <3i 

and the switching times can be found as 

= 1.00268 sec 

t 2 = 2.10893 sec 
t 3 = 3.21518 sec 
tf = 4.21786 sec. 

Note that a longer maneuver time of 4.217 sec is required for a flexible body, as 
compared with the maneuver time of 2.828 sec for a rigidized system. 

The time-optimal control input can then be expressed as 

u(t) = u s (t) - 2u s (t - <j) + 2 u,(t - t 2 ) 

— 2u,(t — t 3 ) + u 3 (t — tf). (8) 
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The open-loop responses of the system to this time-optimal control input are 
shown in Fig. 3 for four different values of k. It can be seen that an ideal, minimum- 
time maneuver is achieved for the nominal system and that perturbations in spring 
constant k, however, result in residual vibrations of the flexible mode. 


3. Preshaped Time-Optimal Control 

As shown in the preceding section, the time-optimal solution requires accurate 
modeling of the system, and the open-loop responses to such an ideal input force 
are very sensitive to the plant modeling uncertainty. In this section, an “impulse- 
sequence” shaping technique developed by Singer and Seering [7-8] is employed to 
preshape the ideal, time-optimal inputs. Performance of the resulting robust feedfor- 
ward controller will be compared with that of a robust, feedback controller in Section 
4. 

The preshaping technique in [7-8] simply utilizes a tapped-delay filter with proper 
weightings and time delays. This technique is briefly reviewed here for our example 
with a single, undamped flexible mode. 

A sequence of m impulses can be expressed in time domain as 

m 

/(0 = X>,<5(f-t.) (9) 

<=i 

with the following normalization 

£>, = i (io) 

1 = 1 

where A, is the magnitude of the i ilx impulse at t = t, and the last impulse occurs at 

t = t m . 

A bang-bang function with (n - 2) switches can be represented as 

«(i) = EW-g (“) 

j=i 

where Bj is the magnitude of a step function at t = t r This bang-bang function ends 
at t = t n . 
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The convolution of u(t ) and /(<) will result in a new multi-switch, multi-level, 
bang-bang function 

n m 

“(0 = - *• ~ *;)• (12) 

•=i j=i 

This function has (mn — 2) switching times and ends at t = ( t m -(- < n ). 

A proper sequence of impulses, whose power spectrum has a notch at a structural 
resonant frequency, can be found as follows. If a sequence of m impulses in Eq. (9) is 
applied to an undamped second-order system with natural frequency of w, the system 
response for t > t m can be represented as: 


m 

E AiU)simj(t — ti ) = Asin{u!t — <f>) 

i=i 

where 


A = 


1 


(5Z A x u>costjt t ) 2 + (^2 AiU>sinu>ti) 2 


i=i 


i=i 


HJT, AiCosuti ' 

If A.’s and s are chosen such that >4 = 0, that is, 


(13) 


A\cosut\ -f Aicosuiti + • ■ • + A m cosut m = 0 (14a) 

A\$inut\ -I- A 2 sinuti + ■ ■ ■ + A m sinu;t m = 0 (14b) 

then the residual vibration will not occur after t = t m . 

Taking derivatives of the preceding two equations for (m — 2) times with respect 
to u>, we get the following 2 (m - 2) robustness constraint equations 

Ai(ti) : sinuti + A 2 {t 2 ) J sinut 2 + • • ■ 

+ A m (t m y sinut m = 0 (15a) 

-4i(fi) J cosu;fi + A 2 (t 2 ) J cosut 2 + • • • 

+ A m {t m y cosut m = 0 (15b) 

j = 1 , • • • , m - 2 . 

For an m-impulse sequence with 1 1 = 0, we now have (2m — 1) equations for (2m — 1) 
unknowns. 
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Figure 4 illustrates three different impulse-sequences with proper A ;' s and the 
time-delay interval of A T = i r/w, where ut is the natural frequency of the flexible 
mode under consideration. 

The frequency response characteristics of this impulse-sequence shaping technique 
can be analyzed simply by taking the Laplace transform of an m-impulse sequence 
as follows: 

L|/W1= 

1=1 

= ( 16 ) 

1=1 

which can be interpreted as a tapped-delay filter (see Fig. 5). The frequency responses 
of this tapped-delay filter for m = 2, 3 and 4 are shown in Fig. 6. It can be seen 
that the frequency component around the resonant frequency is notched out. The 
wider notch width indicates more robustness to frequency uncertainty, but a longer 
response time. 

The flexible-body time-optimal input given by Eq. 8 is now preshaped by a 
tapped-delay filter with m = 2, resulting in the preshaped input command 

u(t) = 0.5u,(f) - u a {t - 1.003) + ti,(< - 2.109) 

+ 0.5u,(< - 2.221) - u s (t - 3.215) 

-«,(< -3.224) + 0.5u,(i -4.218) 

+ u a {t -4.330) -u,(< -5.436) 

+ 0.5u,(t -6.439). (17) 

This preshaped input takes values of ±1.0 and ±0.5. 

The flexible-body time-optimal input is also preshaped using a tapped-delay filter 
with m = 3, resulting in the preshaped control input command 

u(t) = 0.25 n,(0 - 0.5«,(f - 1.003) + 0.5 Uj (* - 2.109) 

+ 0.5u s (f - 2.221) - 0.5 Uj (< - 3.215) 

- u,(t - 3.224) + 0.25uj(< - 4.218) 
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( 18 ) 


+ u,(t - 4.330) ± 0.25u,(f - 4.442) 

- u,(t - 5.436) - 0.5 u a {t - 5.445) 

+ 0.5u,(< — 6.439) ± 0.5u,(< — 6.551) 

- 0.5u,(f - 7.657) ± 0.25u,(< - 8.660). 

It can be seen from Fig. 7 that the preshaped command takes values of ±0.25, ±0.5 
and ±0.75. The time responses of the system to this preshaped input are shown in 
Fig. 8 for four different values of k. It is evident that robustness with respect to 
flexible mode frequency has been increased at the expense of increasing manuever 
time to about 8.66 sec, comparing with the ideal minimum time of 4.218 sec shown 
in Fig. 3. 

The simple, rigid- body time-optimal input given by Eq. 3 can also be preshaped 
to reduce the residual vibration shown in Fig. 2. For example, preshaping the rigid- 

body bang- bang command with a tapped-delay filter with m = 3 results in another 
bang- bang command: 

“(0 = 0.25u,(t) - 0.5u,(t - 1.414) ± 0.5u,(f - 2.221) 

+ 0.25u 4 (< - 2.828) - u,(t - 3.635) 

± 0.25u,(t - 4.442) ± 0.5 u,(t - 5.049) 

- O.Su^f - 5.856) ± 0.25u,(f - 7.270). (19) 

The responses of the system to this preshaped input (not shown here) also indicate 
a reasonable performance robustness (but not better than Fig. 8). In practice, 
however, the time-optimal solution for a rigidized rather than flexible model of a 
multi-link flexible robot (e.g, see [14]) may be preshaped using a tapped-delay filter 
to accommodate the flexible mode effects. 

In summary, we have shown that time-optimal control input, preshaped using a 
tapped delay filter, provides a robust maneuvering scheme which can minimize resid- 
ual structural vibrations. It is also evident that some tradeoffs between performance 
and robustness must always be considered. In the next section, a robust feedback 
compensator will be designed and its performance and robustness will be compared 
with that of the robust, preshaped feedforward approach discussed in this section. 
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4. Robust Feedback Control Design 


As discovered in [10-13], a non-minimum-phase compensation is particularly use- 
ful for practical tradeoffs between performance and robustness for a certain class of 
noncolocated structural control problems. It is, however, often criticized because 
of its sluggish response and its loop gain limitation. In this section, a robust H 0 0 
feedback compensator design is discussed with special emphasis on a proper imple- 
mentation of a non-minimum-phase compensator, incorporating a non-zero set-point 
control scheme. It is shown that a properly designed, feedback controller with a non- 
zero set-point command performs well, as compared with a time-optimal, open-loop 
controller with special preshaping for robustness. 

Consider a single-input single-output (SISO) control system as illustrated in Fig. 9, 
which is the most commonly used configuration for a “two-degree-of- freedom” con- 
troller. The plant and compensator transfer functions are represented as 


I<(s) 

G(s) 


N c (s) 

D c (s) 

Njs) 

D(s) 


(20a) 

(20b) 


where N c (s), D c (s), N(s) and D(s) are polynomials of the Laplace transform variable 
s. The closed-loop transfer function from the desired output command y m to the actual 
output y is then 


yjs) K(s)G(s) 
y-(s) 1 + I<(s)G(s) 1 

NcN 

d c d + n c n f ^ s) ' 


( 21 ) 


Thus, for the conventional feedback control system of Fig. 9, the zeros of the 
closed-loop transfer function are identical with the zeros of the loop transfer function 
K(s)G(s). These zeros sometimes cause an excessive, transient peak overshoot even 
when the closed-loop poles are properly selected. In this case, a prefilter F(s) is often 
used for the cancellation of the undesirable zeros of the closed-loop transfer function. 
(Of course, we cannot cancel the non-minimum-phase zeros!) 
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If the compensator is placed in the feedback path, the closed-loop transfer function 
becomes 


ikl _ ow „ 3l 

Sf(j) 1 + K(s)G(s) ' ’ 

= °' N F(s) 

D,D + NcN ' 1 


( 22 ) 


where the compensator zeros do not appear as zeros of the closed-loop transfer func- 
tion, and a prefilter F(s) must be properly designed for the generation of a control 
input command. Recall that one of the primary motivations for the use of closed-loop 
rather than open-loop control systems in practice is to cope with unexpected distur- 
bances, which an open-loop controller cannot. For this reason, feedforward control is 
seldom used alone, but rather, it is used in combination with feedback control. 

In this section, we show a proper way of implementating a non-minimum-phase 
compensator to minimize such excessive, transient peak overshoot caused by the com- 
pensator zeros. We briefly review a robust H «, control design methodology developed 
in [12-13], and we present a non-zero set-point control scheme for an /f^-based con- 
troller, followed by an example design. 


Robust H a 0 Control 

Consider a linear, time-invariant system described by [15] 


x(t) = Ax(t) -f B\w{t ) + f?2 u (0 

z(t) = C,x(0 + Dnw(t) + Di 2 u(t) (23) 

y(t) — Cix(t) -f D 2 \\v(t) + Z? 22 u (0 

where x(t) is an n-dimensional state vector, w(t) an m^dimensional disturbance 
vector, u(<) an m 2 -dimensional control vector, z(t) a p r dimensional controlled output 
vector, and y(t) a p 2 _< hrnensional measurement vector. 

In order to utilize the concept of an internal feedback loop, the system with 
uncertain parameters is described as [12, 13] 


X 


■ A B x B 2 ' 


X 

z 

— 

C x Dn D12 


w 

. y m 


C 2 D21 D21 


u 
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where C u D u , and D u are not subject to parameter variations. The perturbed 
system matrix in Eq. (24) can be linearly decomposed as follows: 


(25) 


where the first matrix in the right-hand side is the nominal system matrix and Ap is 
the perturbation matrix defined as 


X 

f 

' A 

Bi 

b 2 

z 

= 

c, 

D n 

D , 

1 y . 

1 

c 2 

D 21 

D 2 


1 1 


X 

+ Ap 


w 

J j 


u 


A ^ 

Ap = 


■ AA 

AB X 

ab 2 

0 

0 

0 

<1 
1 

AD 2 i 

AZ?22 


(26) 


Suppose that there are l independent perturbed parameters pi,...,p/ which are 
bounded as |Ap/| < 1. The perturbation matrix A p is then decomposed with respect 
to each parameter variation as 



r m x i 


A P = - 

i 

E 


N , N w JV„ ] = -MEN 


(27) 


where 


E = 


Api 

0 


0 

Ap ( 


(28) 


By introducing the following new variables 

iV r 0 N w N u ] 


2p = 


w p 
XV 

u 


fci n 

Wp — E Zp , 


(29a) 

(29b) 


the perturbed system, Eq. (25), and the input-output decomposition, Eq. (27), can 
be combined as: 


X 


■ A 

M x 

Bi 

B 2 1 

z v 


N x 

0 

N w 

N u 

z 


Ci 

0 

D n 

D 12 

. y . 


. C 2 

My 

D 21 

d 22 . 


x 

w p 

w 

u 


w. 


= — E z. 


(30a) 

(30b) 
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where w p and z v are considered as the fictitious input and output, respectively, due 
to the plant perturbation; and E is considered as a fictitious, internal feedback loop 
gain matrix. 

The following redefinition of z, w, and expansion of the associated matrices enable 
us to employ the state-space representation given by Eq. (23): 


D n 


z J ’ 

Mr B X j , 

0 N w 1 


W < — 



Cx+- 


D\2 *~ 


N x 

Ci 

N u 

D 12 


D 21 


My Z?2i 


( 31 ) 


It can be shown that under certain conditions, there exists an internally stabilizing 
controller such that, for the closed-loop transfer matrix T and for a given design 
variable 7, 

imu < 7 

if and only if the following Riccati equations [15] 


0 = A t X + XA - X{B 2 Bj - \BiBf)X + Ci (32) 

7 

0 = AY + YA t - Y(CjC 2 - ^CjCi)Y + BiBj (33) 

have unique symmetric positive semi-definite solutions A' and Y . 

An //,0-suboptimal controller that satisfies || J^^Hoo < 7, where 7 is a design 
variable specifying an upper bound of the perturbed closed-loop performance T ZX1 . , is 
then obtained as 

u(s) = -h'[sl - Aj^Lyis) (34) 

or 


i c = A c x c + Ly 
- Kx c 


(35a) 

(35b) 


u = 


13 


where 


K = B*X 

(36a) 

H 

1 

* 

” H IT- 

i 

ii 

(36b) 

A c = A + \b x BJX - B 2 I< - LC 2 . 
7 2 

(36c) 


Note that this Hoo controller has a structure similar to a conventional state-space 
controller consisting of an estimator and a regulator, but is designed for a plant system 
matrix 

A+\b x BJX 

T 

where B x B x X/y 2 can be interpreted as an estimate for the worst disturbance input. 
In other words, the separation principle of the conventional linear-quadratic-gaussian 
(LQG) technique does not hold here. Consequently, the non- zero set-point control 
scheme, which has been well established for LQG control synthesis, needs some minor 
modification as discussed in the next section. 

Non-Zero Set-Point Control 

A block diagram representation of a SISO closed-loop system for a conventional 
LQG-type controller with a non-zero set-point command is illustrated in Fig. 10. 
For this configuration, the control input command, u*, corresponding to the desired 
(constant) output command, j/*, is simply given as 

u m = -[C{A - BI\)~ l B]~ 1 y” (37) 

which is independent of the estimator gain matrix L. In this case, we can easily show 
that for dynamic systems having a rigid-body mode, u m depends only on the regulator 
parameters (not on the plant parameters such as m x , m 2 , and k of the example model 
shown in Fig. 1). Hence, the non- zero set-point control scheme is inherently robust 
to plant parameter uncertainty for a certain class of dynamical systems with at least 
one pole at the origin (i.e., a type-1 system). 

Now consider a closed-loop control system with an Hqq controller as shown in Fig. 
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11, which is described in state-space form as 
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A -BiK 1 
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’ b 2 ' 

. ic . 
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+ 
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C 2 0 
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X c 


(38a) 


The input command u* corresponding to the desired output y* can be simply found 
as 

u* = -[CA-'B ]- V 


where 


A = 

1 1 

-b 2 k ' 

A c 

, B = 

' b 2 
. b 2 

C = 

C 2 0 

1- 




(39) 


(40) 


Similarly to the LQG case, it can be shown that for dynamic systems having a 
rigid-body mode, u* for an H ^ controller depends only on the controller parameters 
(not on the plant parameters such as m lt m 2 , and k of the example model shown in 
Fig. 1). However, u* now depends on both the gain matrices K and L , not just on 
the regulator gain matrix K as for the LQG case. 


Example Design 

We now consider the two-mass-spring model shown in Fig. 1. A control input 
force u acts on body 1 and the position of body 2 is measured as y, resulting in 
the so-called noncolocated control problem. The control design objective here is to 
achieve a fast settling time (high performance) for an output command y' = 1 and 
robust performance over a range of spring stiffness uncertainty considered in Section 
3. The control input is bounded as |u| < 1 and the system has the nominal values of 
mj = m 2 = k = 1. 

The plant model can be represented in state-space form as: 
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(41a) 
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y = x 2 + v 
z = x 2 


(41b) 

(41c) 


where u?i and v are, respectively, plant disturbance and sensor noise assumed for 
control design purposes. 

The variation AA is decomposed as 


A A = -A k 


’ 0 
0 
1 

-1 


[ 1 -1 0 0 . 


(42) 


The other elements of the perturbation matrix in Eq. (26) are all zeros. Note that 
AA is spanned by the matrices 

’ 0 

Mx= J 
. -1 

where M x is the fictitious disturbance distribution matrix spanning the columns of 
Aj 4, and N x is the fictitious controlled output distribution matrix spanning the rows 
of AA The fictitious input and output for this example are expressed as 


and N x = [ 1 -1 0 0 ] 


Zp = Nx = Xi — x 2 and w p = —A k z v . 


(43) 


Equation (43) replaces the parameter variation in Eqs. (41), resulting in the fol- 
lowing equations: 
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(44b) 
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y = x 2 + V 


(44c) 


where k is the nominal spring constant. As can be seen in Eq. (44b), control input u 
is also included in the performance variable z in order to minimize control effort. 

By defining 


w 


w v 

Wi 

V 


and by using the definition of Eq. 
represented as 


and z 


(45) 


31), the system matrices in Eq. (23) can be 


0 0 10 ' 
0 0 0 1 
-k k 0 0 
k -k 0 0 
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1 1 0 

1 


1 

1 

\ — ‘ 
0 

0 

1 


. 0 . 



■ 1 

-1 0 

0 ' 


' 0 ' 

<?! = 

0 

1 

0 

0 
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_ 0 
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0 
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1 

c 2 = | 

0 

1 0 0 

]. 

D 21 = 

[» 


and D\\ — 03*3, D 22 — 0. 

For the example design considered here, the disturbances w p , w\, and v are mul- 
tiplied by weighting factors, 0.1, 0.025, and 0.025, respectively. The performance 
specification bound 7 is chosen to be 1. The weighting factors and 7 represent rela- 
tive disturbance levels and overall closed-loop performance level, respectively. 

By solving two Riccati equations, Eqs. (32) and (33), we get the controller gain 
matrices K and L as follows: 


K = 
L = 


1.506 -0.494 1.738 0.932 ] , 
0.720 2.973 -3.370 4.419 ] T . 
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The corresponding robust Hoo controller is then 

0.0827(5^5 + lK-ak 


u(s) = - 


+ 1 ) 


[(rfe) ! + 2 (0-825)(rfe) + i! 
(ifc + i) 


Krfs ) 2 + 2(0.459)( 5 * 5 ) + 1] 


ttivM 


( 47 ) 


which is a non-rninimum-phase compensator. 

Figure 12 shows a closed-loop root locus versus overall gain of this compensator. 
The gain margin is 3.28 dB, and the closed-loop system is stable for 0.44 < k < 3.27, 
which corresponds to -0.56 < A k< +2.27. The nominal system’s closed-loop poles 

are: 

-0.337 ±0.336j, -0.514 ± 0.414j, (48) 

-0.376 ±1.495;, -1.109 ± 1.797;. 

The time response of the closed-loop system implemented as in Fig. 9, for F(s) — 1 
and y m = 1, is shown in Fig. 13. A non-minimum-phase behavior of the closed-loop 
system is evident and the nominal system has a peak overshoot of about 80% and a 
settling time of 15 sec. When compared with the response of a feedforward controller, 
shown in Fig. 8 , the overall response is not acceptable. It can be shown that the 
excessive overshoot is due to the compensator zero at s = —0.145. This zero may 
be cancelled by a prefilter F{s), but the resulting slower settling time may not be 
desirable. 

Figure 14 shows the closed- loop responses for four different values of k to an 
output command of y’ = 1 (consequently, u* = 0.9959), when the same controller is 
implemented as in Fig. 11. Clearly, the responses no longer have excessive overshoot 
and the settling time is quite short, as compared with the response in Fig. 13. The 
overall responses are also comparable with those of a preshaped feedforward controller, 
as can be seen in Fig. 14. The control input u{t) is always within the saturation limit 

of one. 

In summary, we may conclude that a feedback controller, when implemented prop- 
erly, could achieve good performance and robustness, for both command following as 
well as disturbance rejection problems. The proposed feedforward/feedback control 
approach is robust for a certain class of uncertain dynamical systems, since the con- 
trol input command computed for a given desired output does not depend on the 
plant parameters. 
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5. Conclusions 


In this paper we have investigated both feedforward and feedback control ap- 
proaches for rapid maneuvering control of uncertain dynamical systems. A simple 
two-mass-spring example was used to illustrate the control concepts and methodolo- 
gies. It was shown that a time-optimal control input, preshaped using a tapped-delay 
filter, provides a rapid maneuver and robust suppression of residual structural vibra- 
tions. A proper implementation of a non-minimum-phase compensator with a non- 
zero set-point control command was discussed. It was demonstrated that a properly 
implemented feedback controller performs well, when compared with a time-optimal, 
open-loop controller with special preshaping for performance robustness. 
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Figure 2: Responses to rigid-body, time-optimal control input. 
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Figure 3: Responses to flexible-body, time-optimal control input. 
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Figure 4: Impulse-sequence shaping. 




Figure 5: Tapped-delay filter. 












Figure 8: Responses to preshaped time-optimal command# 
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Figure 9: Conventional feedforward/feedback control system configuration. 







Figure 10: Control system configuration with LQG-type controller. 








Figure 11: Control system configuration with //<»- type controller. 






Figure 12: Root locus versus overall loop gain. 
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Figure 13: Time responses for conventional control configuration 
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Abstract 

A new approach for computing time-optimal open-loop control inputs for uncertain 
flexible spacecraft is developed. In particular, the single-axis, rest-to-rest maneuvering 
problem of flexible spacecraft in the presence of uncertainty in model parameters is 
investigated. Robust time-optimal control inputs are obtained by solving a parameter 
optimization problem subject to robustness constraints. A simple dynamical system 
with a rigid-body mode and one flexible mode is used to illustrate the concept. 


1. Introduction 

This paper is concerned with the problem of computing open-loop control inputs for 
flexible spacecraft, robotic manipulators and pointing systems in space, which are often 
required to maneuver as quickly as possible without significant structural vibrations 
during and/or after a maneuver. 

A standard, time-optimal control approach to such a problem requires an accurate 
mathematical model, and thus the resulting solution is often sensitive to variations in 
model parameters. For this reason, an open-loop time-optimal controller is seldom used 
in practice. Consequently, the development of a “robustified” open-loop approach for a 
rapid maneuver without significant structural vibrations is of current research interest 
[1-3]. Other open-loop approaches [4-6] attempt to find a smooth continuous forcing 
function (e.g., a versine function) that begins and ends with zero slope. The basic idea 
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tGraduate Research Assistant, Student Member AIAA. 
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behind such approaches is that a smooth control input without sharp transitions is less 
likely to excite structural modes during maneuvers. 

In this paper, a new approach is developed for computing time-optimal control inputs 
for the single-axis, rest-to-rest maneuvering problem of flexible spacecraft in the presence 
of structural mode frequency uncertainty. A parameter optimization problem, where 
the objective function to be minimized is the maneuvering time, is formulated with 
additional constraints for robustness with respect to structural parameter uncertainty. 
The resulting robust time-optimal solution is a multi-switch bang-bang control which 
can be implemented for spacecraft equipped with on-off reaction jets [7], This result 
further confirms that most open-loop approaches, which utilize a smooth continuous 
control input so that structural modes are less likely to be excited, do not fully utilize 
the available control energy in performing a robust time-optimal maneuver. 

The paper is organized as follows. In Section 2, the rest-to-rest maneuvering con- 
straints for multi-switch bang-bang inputs are derived with some discussion on the pre- 
vious results of [8-10] on the number of switchings for the time-optimal solution. The 
standard, time-optimal control problem is then transformed into a constrained parameter 
optimization problem. In Section 3, robustness constraints are derived and incorporated 
with the parameter optimization problem formulated in Section 2. A simple dynamical 
system with a rigid-body mode and one flexible mode, shown in Fig. 1, is used to illus- 
trate the concept, and the robust time-optimal solution for a case with a single two-sided 
control (Case 1) is discussed. In Section 4, the same example with two one-sided control 
inputs (Case 2) is further investigated. 

2. Time-Optimal Rest-to-Rest Maneuver 
Problem Formulation 

Consider a linear model of flexible spacecraft described by 

Mx + Kx = Gu (1) 

where x is a generalized displacement vector, M a mass matrix, K a stiffness matrix, G 
a control input distribution matrix and u a control input vector. 

In this section, we consider a case with a scalar control input u(t) bounded as 

-1<u<1 (2) 
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Equation (1) can be transformed into the decoupled modal equations: 

Vi +ufyi = fau 
y 2 + w|j /2 = 4>7 u 

Vn + = 4> n u 


where y { is the i th modal coordinate, u \ the i th modal frequency, fa i tk modal gain, and 
n the number of modes considered in control design. 

The problem is to find the control input which minimizes the performance index 


subject to Eqs. (2) and (3), and given boundary conditions. 

The time-optimal control problem of a linear controllable system has a unique solu- 
tion which is “bang- bang” control with a finite number of switches [10]. For a spring-mass 
dynamical system with n degrees of freedom, the time-optimal bang-bang solution for 
a rest-to-rest maneuver has, in most cases, (2n — 1) switches [8-9], and the solution is 
symmetric about tj/ 2. That is, for a case with (2n — 1) switches, we have 


— * ^2n ^2n — J l?"'"i^ 


(4) 


where t 2n = tj- 

A bang-bang input with (2n - 1) switches can then be represented as: 

u(t) = ZBj«.(t-tj) W 

J=0 

where B : is the magnitude of a unit step function u 3 (t) at tj. This function can be 
characterized by its switch pattern as: 

B = { Bo, B \ , B2, • • * B%n } 

T = { U, i'll * * * y ^2 n } 

where B represents a set of B : with Bq = B 2 n = il and Bj = i2 for j = 1, • • • ,2n — 1, 
T represents a set of switching times (t\, • * • ,^ 2 n-i) an< l the initial and final times (to = 0 

and tf = t 2 n)< 
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Rest-to-Rest Maneuver Constraints 


Consider the rigid-body mode equation with uj\ = 0: 

Vi = 4>\u ( 6 ) 

with the rest-to-rest maneuvering boundary conditions 

Jh(0) = 0, yi{t/) ^ 0 

3/i (0) = 0, y 1 (t / ) = 0 (7) 

Substituting Eq. (5) into Eq. (6) and solving for the time response of the rigid-body 
mode, we get 

yx(t > t}) = y Y;(tf -ttfB, (8) 

j=0 

The rest-to-rest maneuvering constraint for the rigid-body mode can then be written as: 


4 £(1/ - -»,(</) = 0 (9) 

J=° 

Consider now the structural modes described by 

Vi + u*y, = i = 2, -,n (10) 

with the corresponding boundary conditions for the rest-to-rest maneuver: 

y.(0) = 0, y,(tj) = 0 

y.(0) = o, yi(tf) = o (I1) 

for each flexible mode. 

Substituting Eq. (5) into the i th structural mode equation and solving for the time 
response for t > tf, we get 


d> ^ n 

y*(0 = — j^2^j cosuj .( t -tj) 

u ' J=0 

S' 2n 

= j[cosu,(t cos -t n ) 

U i j=0 

2 n 

+ siniL>,(t - t n )Y, - f n )] 

j= o 

It can be shown that the following constraint equation for each mode 


2 n 


^BjSinu;,^ - <„) = 0 
]=0 


( 12 ) 


(13) 
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is always satisfied for any bang-bang input which is symmetric about the mid-maneuver 
time t n . Consequently, we have the following flexible mode constraints for no residual 
structural vibration (i.e., y,(t) = 0 for t > tj ): 

2 n 

^ B } cos u,(tj - t n ) = 0 (14) 

j=0 

for each flexible mode. 

Parameter Optimization Problem 

For a spring-mass system of n degrees of freedom, the time-optimal solution repre- 
sented by Eq. (5) has the (2n - 1) unknown switching times and the final time t f to be 
determined. The time-optimal control problem can now be formulated as a constrained 
parameter optimization problem as follows: 

Determine a control of the form given by Eq. (5) that minimizes the final time tj 
subject to 


T Blz-'iM -».('/) = ° 

z J=0 

2n 

Y', Bj cosu >j{tj — t n ) = 0; i = 2, •••,n 

j=0 

tj > 0; j = l,---,2n 

where tj = t^n- 

Remark: Note that <f>i (i = 2, • • • , n) do not appear in Eqs. (15); that is, the optimal 
solution to this problem is independent of the flexible mode shapes. In other words, the 
time-optimal control input is independent of actuator location for a system described 
by Eq. (3) with a scalar input. 

Standard optimization packages (e.g., IMSL subroutines) can be used to obtain the 
solution of the above optimization problem. The major advantage of the proposed 
approach, compared to other direct numerical optimization approaches employed in 
[11-14] for the time-optimal control problem, is that some robustness constraints with 
respect to plant parameter uncertainty can be easily augmented. This subject will be 
discussed in detail in Section 3. 


(15a) 

(15b) 
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Sufficient Conditions for Optimality 


Equations (15) are necessary conditions for time-optimal control, and sufficient con- 
ditions for optimality can be checked as follows. 

Let the costate vector corresponding to the modal state vector [t/! , y, , uy 2 ,y 2 , • • • , u n y n y n ) T 

be defined as 

)T 


M0=[M0» ?i(0. Pn(t), q n (t) j J 

It is shown in [8] that at mid-maneuver, we have 


(16) 


^(fn) — [ Pl(fn)t 0, Piitn), 0, •••, p n (t n ), 0 J 

and thus the costate vector can be solved as 

Pi(0 = Pi(*») 

?i(0 = -(< ~ U)P\{t n ) 

Pi(0 = Pi(t„) cos u;,(t - t„ ) 

9.(0 = -Pi(f n )sinw,-(f -t„) 

for each flexible mode. Then p,{t n ), for i = 1, • • • , n, can be found from the following n 
linear equations: 

n 

1 + 4>iPi(t n )t n + £ <t>iPx{t n ) sin Uit n = 0 (19) 

i=2 

n 

0lPl(fn)(<i - tn) + £4>,p,(t n )silUJi(tj -t„) = 0 (20) 

1 1=2 
where j=n + l,---,2n-l. 

The solution obtained by minimizing tj subject to Eqs. (15) becomes time-optimal 
provided that 


(17) 


(18) 


S(t) - -4>\P\(t n )(t - t n ) - <£,?,(*„) sin u>,(< -i n ) ^ 0 (21) 

i=2 

for t 6 (t n , t 2n ) and t ^ tj, j = n + 1, • ■ ■ , 2n, and S(t) represents the switching function. 

Remark: If modal frequencies are rational multiples of each other and the funda- 
mental frequency, u> 2 , satisfies the following relationship 


u 2 = 2£ir 




4>x 


y\(h) 


£ = 1,2, 


( 22 ) 


then the time-optimal solution has only one switch and is equivalent to the solution of 
a “rigidized” case. 
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Example: Case 1 with a Scalar Control Input 

Consider a simple example, shown in Fig. 1, which is a generic representation of a 
flexible spacecraft with a rigid-body mode and one flexible mode. Case 1 with a scalar 
control input u(t ) is considered here. The equations of motion are 

miii + k(x i — x 2 ) = ti t = u (23a) 

m 2 x 2 — k(ii — ^ 2 ) = u 2 = 0 (23b) 


where x, and x 2 are the positions of body 1 and body 2, respectively, and the nominal 
parameters are m x = m 2 = k = 1 with appropriate units, and time is in units of second. 
The boundary conditions for a rest-to-rest maneuver are given as 

xi(0) = i 2 (0) = 0, xi(tj) = x 2 (tj) = 1 ( 24 ) 

x,(0) = x 2 (0) = 0, i,(t/) = i 2 {tj) = 0 

The modal equations are 

y"i = u/2 (25a) 

y 2 + u 2 y 2 = u/2 (25b) 

where u = \/2 rad/sec is the nominal flexible mode frequency. The corresponding 
boundary conditions for modal coordinates are 

2 /i( 0 ) = 3 / 2 ( 0 ) = 0 , y\{tf) = 1 , i/ 2 (*/) = 0 (26) 

3 /i( 0 ) = 3 / 2 ( 0 ) = 0 , y\(tj) = Viitf) = 0 

Since there are 3 switches, the time-optimal switch pattern for the given boundary 
conditions is represented as 

B={B 0 , B u B 2 , B 3 , B<} 

= { 1 , - 2 , 2 , - 2 , 1 } 

T = { to, t\, t 2 , < 3 , <4 } 

with the symmetry conditions 


U = 2< 2 

Z3 = 2 1 2 — 1 1 
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The time-optimal control problem is then formulated as the following constrained 
minimization problem: 


subject to 


min J = 2*2 


(27) 


2 + 2*j + *2 — 4*i*2 = 0 

1 — 2 cosu ;(*2 — *i) + cosu ;*2 = 0 

*i, *2 > 0; 

A standard IMSL FORTRAN subroutine was used to obtain a solution as: t x = 1.003 
and < 2 = 2.109. The computed solution satisfies the optimality conditions of Eq. (21); 
i.e., the switching function vanishes only at * = t u t 2 and * 3 . Thus, the solution obtained 
via Eq. (27) is indeed time-optimal and it can be expressed as 

«(<) = «.(<) - 2 u,(* - 1.003) + 2 u 3 (t - 2.109) 

— 2u # (* — 3.215) + u 3 (t — 4.218) (28) 

The time responses of X 2 to the time-optimal control input are shown in Fig. 2 for 
four different values of k. It can be seen that the resulting responses are sensitive to 
variations in the model parameter k. 

3. Robust Time-Optimal Control 

As discussed in the preceding section, a standard, time-optimal control approach 
requires an accurate mathematical model and thus the resulting solution is often sensitive 
to plant modeling uncertainty. 

In this section, a new approach, expanding on the approach introduced in Section 2, 
is developed for computing time-optimal control inputs for the single-axis, rest-to-rest 
maneuvering problem of flexible spacecraft in the presence of structural mode frequency 
uncertainty. A parameter optimization problem, where the objective function to be 
minimized is the maneuvering time, is formulated with additional constraints for ro- 
bustness with respect to the structural frequency uncertainty. The resulting robustified, 
time-optimal soultion is a multi-switch bang-bang control, and is thus implementable 
for spacecraft equipped with on-off reaction jets [7]. 
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Robustness Constraints 


By taking the derivative of Eq. (12) with respect to we get 


<Lo x 


^icosw,(< - y)E(^ - - *-+) 


for each flexible mode. Letting dy^tj/tLi, = 0 for all t > tj , we have 


( 29 ) 


f(l j -|)B jS inu,,(( J -^) = 0; i = 2, • ■ • , n 


(30) 


which is called the first-order robustness constraints. 

Similarly, taking the derivative of Eq. (12) r, times with respect to u>, results in r‘ h 
order robustness constraints for each flexible mode as follows: 


D< i-^) m B J sin,v,(i ) -^) = 0 
)- 0 1 1 

for m = 1,3, rj (31a) 

cosat.ll, - |) = 0 

J=0 2 

/or m = 2, 4, • • • < r, (31b) 


There are total r robustness constraints for (n — 1) flexible modes, where 

r = ±r t (32) 

i=2 

If these robustness constraints are included in the constrained minimization problem 
formulation described by Eq. (15), the number of switches in the bang-bang control 
input, in most cases, must be increased to match the number of the constraint equa- 
tions. Due to the symmetric nature of the rest-to-rest maneuvering problem, adding one 
robustness constraint will require, at least, two more switches. 

Robust Time-Optimal Control 

If r robustness constraints are considered for a flexible spacecraft of n modes, the 
corresponding robustified bang-bang control input becomes 

2(n+r) 

«(()= E B , ».(‘-‘i) < 33 > 

j=0 
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which has 2 (n + r) unknown switching times. Due to the symmetry property of the 
optimal solution for the rest-to-rest maneuvering problem, we have 

i i = *2(n+r) - <2(n+r)-j; j = 1, * * • , M + r (34) 

Therefore, there are only (n+r) unknowns to be determined in Eq. (33). These unknowns 
can be determined by minimizing t f subject to the (n + r) constraint equations: one 
positioning constraint for the rigid-body mode, (n - 1) no-vibration constraints and r 
robustness constraints. 

While many theoretical issues (e.g., the uniqueness of the optimal solution) need to 
be explored, a solution can be obtained by solving the following constrained parameter 
optimization problem: 

min J = tj = < 2 (n+r) (35) 

subject to 

<t>\ 

T ^1 (*2(n+r) — t : ) 7 Bj — Vl(tf) = 0 

6 1=0 
2(n+r) 

XI Bj cos Uiitj - t n+T ) = o 
1=0 
2(n+r) 

XI (^7 — ^n+r) Bj sinu?j(<y — t n+r ) = Oj 

J=0 

for m = 1,3, ••• < r, 

2(n+r) 

XI (^/ — ^n+f) Bj COSU li(tj — t„+ r ) = 0; 

1=0 

for m = 2,4, • • ■ < r, 

> 0; j = 1, ■ ■ • ,2(n + r) 

for each flexible mode. The resulting bang-bang control input, which has 2r more 
switches than the time-optimal bang-bang solution of Section 2, is called a robust (or 
robustified) time-optimal solution in this paper. 

Example: Case 1 with a Scalar Control Input 

For Case 1, the time-optimal control is a three-switch bang-bang function, but the 
resulting responses were shown to be sensitive to variations in model parameter k . A 
robust time-optimal solution of the same problem is now developed as follows. The 
switching pattern for a case with the first-order robustness constraint is assumed as: 
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B = { Boy B \ , Bs, i?*, i?5, B$ } 

= { 1, -2, 2, — 2, 2, -2, 1} 

(36a) 

T = { toy t\y ^3) t'A'i ^5? ^6 } 

(36b) 

with the symmetry conditions 

< 4 = 2<3 — 

< 5 = 2<3 — t\ 

(37) 

— 2<3 

The constrained optimization problem with the first-order robustness constraint can 

be formulated as: 

min J = t 6 

(38) 

subject to 

2 + 2 1]- 2 1\ -t\- 4<i< 3 + 4t 2 < 3 = 0 

(39a) 

COSW<3 — 2cOSO>(t3 — <i) 

+ 2 cosu;(f3 — t 2 ) — 1=0 

(39b) 

< 3 sinu>< 3 - 2(<3 - U) sinw(< 3 - h ) 
+ 2(<3 - < 2 ) sinu )(t 3 - t 2 ) = 0 

(39c) 

*i) h, h > 0; 

A robust time-optimal solution with 5 switches can be found as: 

h = 0.7124, < 2 = 1.6563 
< 3 = 2.9330, t 4 = 4.2097 

(40) 

<s = 5.1536, t 6 = 5.8660 



The time responses of x 2 to this “robustified” time-optimal control input are shown 
in Fig. 3 for four different values of k. It can be seen that the resulting responses are less 
sensitive to parameter variations, compared to the responses to the ideal, time-optimal 
control input as shown in Fig. 2. Performance robustness has been increased at the 
expense of the increased maneuvering time of 5.866 sec, comparing to the ideal minimum- 
time of 4.218 sec. It is, however, emphasized that simply prolonging the maneuver time 
does not help to reduce residual structural vibrations caused by modeling uncertainty. 
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Remarks 

An impulse-sequence shaping technique, developed by Singer and Seering [1,2], was 
employed by Wie and Liu [3] to preshape the ideal, time-optimal control input given in 
Eq. (28). For example, the two-impulse preshaped bang-bang command was obtained 

as 


u(t) = 0.5u,(t) — u,(t — 1.003) + u 3 (t — 2.109) 

+ 0.5u a (<- 2.221) 3.215) 

- u 3 {t - 3.224) + 0.5u,(< - 4.218) 

+ u 3 (t - 4.330) - u 3 (t - 5.436) 

+ 0.5u,(t — 6.439). ( 41 ) 

This preshaped input takes values of ±1.0 and ±0.5, and the resulting response becomes 
less sensitive to flexible mode frequency variations as demonstrated in [3]. 

For Case 1, the time responses to the time-optimal input of Eq. (28), the robust 
time-optimal input of Fig. 3, and the preshaped inputs in [3] can be compared as shown 
in Figs. 4 and 5. Comparisons of the maneuvering times of four different control schemes 
are illustrated in Fig. 4. Parameter robustness with respect to spring constant variations 
is compared in Fig. 5 . From these figures, it is evident that the proposed approach of this 
paper provides a faster and more robust maneuver than other robustified feedforward 
approaches. Also, unlike other approaches in [1-6], the resulting solution of our new 
approach is a multi-switch bang-bang control which can be implemented for spacecraft 
with on-off reaction jets. 

In the next section, we will consider a case with two one-sided control inputs in order 
to further explore a u time-optimal” actuator placement problem. 

4. A Case with Two One-Sided Control Inputs 
Problem Formulation 

Consider Case 2, illustrated in Fig. 1, with two one-sided control inputs bounded as 

0 < ui < +1 ( 42a ) 

- 1 < u 2 < 0 (42b) 
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Since the control inputs are one-sided, each control input for the time-optimal solution 
need not be an odd function about the mid-maneuver time. Thus, the problem with 
one-sided control inputs becomes more difficult to solve than the standard problem with 
two-sided control inputs, and many theoretical issues (e.g., the uniqueness and structure 
of time-optimal solutions) need further investigation. 

For Case 2, the modal equations of the system with nominal parameter values are 


Vi — ^ Ul + u 2) 
2 1, 

y 2 + w y 2 = -(u, - u 2 ) 


(43a) 

(43b) 


where u> = \/2 rad/sec is the nominal flexible mode frequency. 

For the control input constraint given by Eq. (42), the control inputs can be expressed 
as 


/v-i 

U 1 = E M* -tj) -u 3 (t -tj - Aj)] (44a) 

0,2,4,- 

N 

u 2 =- £ lu s (t -tj) -u a (t -tj - Aj)] (44b) 

i=1.3,5,- 

which is in the form of one-sided pulse sequences as shown in Fig. 6. The j ifl pulse starts 
at tj and ends at (tj + Aj). Due to the symmetric nature of the rest-to-rest maneuvering 
problem, we assume that u t and u 2 have the same number of pulses, (N + l)/2, where 
N is defined as in Fig. 6. 

Substituting Eq. (44) into Eq. (43a) and solving for the time response of the rigid- 
body mode, we get 


»■(< S */) - X D-n'plA, - 2i,A y - AJ] (45) 

* >=o 

For the desired boundary condition, yi(< > </) = 1, the following constraint must hold 


£(— 1)'A, = 0 (46) 

j=° 

The positioning constraint for the rigid-body mode then becomes 

B-1)'[2(,A j +A;|+4 = 0 (47) 

j=0 
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Substituting Eq. (44) into Eq. (43b) and solving for the time response of the flexible 
mode, we get 

1 N 

y 2 (t) = — cos ut [cos utj — cosu(tj + Aj)] 

4 ,=o 

1 N 

smut V'[sinu;f, — sinu;(*j + Aj)] (48) 

4 i-o 

for t >tj. 

Also, rest-to-rest maneuvering requires y 2 (t) = 0 for t > t/; i.e., we have 

N 


y)fcos utj — cos u{tj + Aj)] = 0 
1=0 

(49a) 

J — ^ 

N 

^[sinwij — sinu;(tj + A ; )] = 0 

(49b) 


3=0 


which become the no-vibration constraints for the rest-to-rest maneuvering problem. 

Time-Optimal Control 

Let the time-optimal control inputs for the rest-to-rest maneuver problem be of the 
form 

tii = «»(<) - t i,(t ~ A) 

u 2 = — u,(t — <i) + u,(t — t 2 — A) 

where each input has a single pulse with the same pulse width of A, is defined as 
shown in Fig. 6, and the maneuver time tj = t\ 4- A. 

The rest-to-rest maneuver constraints can be obtained from Eqs. (47) and (49) as: 

tj — (2/A) - A = 0 (50a) 

sin(u;t//2) + sin(u>(A — tj/ 2)) = 0 (50b) 

which can be combined as: 

sin(u>A/2) cos(u;/(2A)) = 0 (51) 

The time-optimal solution can then be obtained by solving the constrained mini- 
mization problem: 

min J = tf = (2/ A) -F A (52) 
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subject to the constraint given by Eq. (51). 

The solution of this problem can be found as 

A = 0.9003 
t f = 3.1218 

The time responses of x 2 to the time-optimal control inputs are shown in Fig. 7 
for four different values of k. The maneuver time and control on-time are respectively 
3.12 sec and 1.8 sec. As expected, the resulting responses are sensitive to parameter 
variations. 

Remark: A most interesting feature of this solution is that the overall input shape 
shown in Fig. 7 is of a “bang-off-bang” type, resulting in the control on-time of 1.8 sec 
which is different from the maneuver time t f of 3.12 sec. For Case 1 and a rigidized case 
[3], both the maneuver time and control on-time are 4.218 sec and 2.828 sec, respectively. 
Therefore, the actuator configuration for Case 2 is considered to be “optimal” in the 
sense of minimizing both the maneuver time and control on-time. 

Robust Time-Optimal Control 

Similar to Case 1 in Section 2, we now consider the “robustification” of the time- 
optimal solution obtained in the preceding section. 

Letting the derivative of Eq. (48) with respect to u> be zero, we get 

dy 2 1 

— = - -sinwf cosujfj - (tj + AjJcosu^t, + A ; )] 

UUI * J=0 

1 N 

+ - cos ut sin utj -(tj- |- Aj) sinw^ + A,)] 

J=0 

= 0 (53) 

For this derivative to be zero for arbitary t > t f , we must have 

N 

cosu>t] — (tj + Aj) cosu }(tj + Aj)] = 0 (54a) 

j = o 

N 

sin ut } - (tj -I- A J -)sinw(t j + A,-)] = 0 (54b) 

i=o 

which is called the first-order robustness constraints. 
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Taking the derivative of Eq. (48) r times with respect to u>, we get the r^-order 
robustness constraint equations for input pulse sequences as follows: 

5Z[(*j) m cos u)t } — (tj + Aj) m cos u>(tj + Aj)] = 0 (55) 

j=o 

N 

^[(<j) m sinu;<j — (tj + Aj) m sinu;(<j + Aj)] = 0 (56) 

:= o 

for m = 1,2, • • • , r. 

As an example, we consider the first-order robustness constraint, incorporated with 
the rest-to-rest maneuver constraints, to construct robust time-optimal pulse sequences. 
Assuming that each input has two pulses, we can represent the control inputs as follows. 

tti = u a (t) — u s (t — Ao) + u a (t — t 2 ) 

— u a (t — t 2 — A 2 ) (57a) 

u 2 = — u,(t — 1\) + u,(t — t x — Ai) — u a (t — < 3 ) 

+ u a (t — t 3 — A 3 ) (57b) 

in which we have seven unknowns to be determined, and tj and Aj are defined as shown 
in Fig. 6. 

The robust time-optimal solution can then be obtained by solving the constrained 
parameter optimization problem 

min J = t 3 + A 3 (58) 


subject to 

Ao — Ai + A 2 — A 3 4- A 4 = 0 
£(-1)'[2( j A, + A;] + 4 = 0 

3 

^[coscj/j — cos uj(tj + Aj)] = 0 

j=0 

3 

^[sina;<j — sin u>(tj + Aj)] = 0 
;=o 
3 

^[tjCOSLJtj - (tj + Aj)cosu;(<j + A ; )] = 0 
i=o 
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3 

£(*•> sinu;< ; - {t } + A^sinw^ + A,)] = 0 

3=0 

Aj>0; > = 0, 1, 2, 3 

U, hi h > 0 

The solution to this problem can be obtained as: 

<o = 0.0000, A 0 = 0.4274 
= 2.3357, A, = 0.4329 

h = 2.1132, A 2 = 0.4329 ( 60 ) 

< 3 = 4.4544, A 3 = 0.4274 

The time responses of x 2 to the robust time-optimal control inputs are shown in 
Fig. 8 for four different values of k. It is seen that the robustness has been increased at 
the expense of the increased maneuvering time of 4.882 sec, comparing with the ideal 
minimum-time of 3.122 sec. However, note that the control on-time is only 1.721 sec, 
compared to the control on-time of 1.8 sec of the ideal, time-optimal solution. 

6. Conclusions 

A new approach to the robust time-optimal control of uncertain flexible spacecraft 
has been investigated. The unique feature of the proposed approach is the fairly straight- 
forward incorporation of the robustness constraints into a standard parameter optimiza- 
tion problem, where the objective function to be minimized is the maneuvering time. 
The case with two one-sided control inputs has shown an interesting feature from the 

viewpoint of “optimal” actuator placement for minimizing both the maneuver time and 
control on-time. 
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Case 1: -1<u,<1; u 2 =0 
Case 2: 0 ^ u ! ^ 1 ; -1 ^ u 2 - 0 


Figure 1: Generic model with a rigid-body mode and one flexible mode. 
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Figure 2: Responses to time-optimal control input (Case 1). 
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Figure 3: Responses to robust time-optimal control input (Case 1). 
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Maneuver Time t f 



Figure 4: Comparision of maneuvering time. 
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Residual Vibration Amplitude for y i (t f ) = i 


1- Time-optimal bang-bang (Fig. 2) 

2- Robust five-switch bang-bang (Fig. 3) 

3- Two-impulse preshaped bang-bang (Ref. 3) 

4- Three-impulse preshaped bang-bang (Ref. 3) 


Spring Stiffness k 


Figure 5: Comparision of parameter sensitivity 
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x2 



Figure 7: Responses to time-optimal control inputs (Case 2). 
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Figure 8: Responses to robust time-optimal control inputs (Case 2). 
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Abstract 

The problem of computing open-loop, on-off jet firing logic for flexible spacecraft in 
the face of plant modeling uncertainty is investigated. The primary control objective 
is to achieve a fast maneuvering time with a minimum of structural vibrations during 
and/or after a maneuver. This paper is also concerned with the problem of selecting a 
proper pair of jets for practical trade-offs among the maneuvering time, fuel consump- 
tion, structural mode excitation, and performance robustness. A time-optimal control 
problem subject to parameter robustness constraints is formulated. A three-mass-spring 
model of flexible spacecraft with a rigid-body mode and two flexible modes is used to 
illustrate the concept. 
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1. Introduction 


This paper deals with the problem of computing open- loop, on-off jet firing logic for 
flexible spacecraft which are sometimes required to maneuver as quickly as possible with 
a minimum of structural vibrations during and/or after a maneuver. Most standard time- 
optimal control approaches to such a problem require an accurate mathematical model, 
and thus the resulting solution becomes sensitive to variations in model parameters. 

Expanding on the recent results of [1-2], we further explore the robust time-optimal 
control problem of flexible spacecraft in the face of modeling uncertainty. In particular, 
we study the problem of selecting a proper pair of jets for practical trade-offs among 
the maneuvering time, fuel consumption, structural mode excitation, and performance 
robustness. A parameter optimization approach, with additional constraints for per- 
formance robustness with respect to modeling uncertainty, is employed to solve such 
a robust time-optimal control problem. However, many theoretical and practical im- 
plementation issues inherent to constrained parameter optimization problems are not 
elaborated in this paper. 

A simple math model of flexible spacecraft with a rigid-body mode and two flexible 
modes, as shown in Fig. 1, is used to illustrate the concept and methodology. We consider 
the case in which the structural flexibility and mass distribution of the vehicle are quite 
uncertain, while the total mass (or inertia) of the vehicle is well known. Consequently, we 
focus on the robust control problem of flexible spacecraft in the face of modal frequency 
uncertainty as well as mode shape uncertainty. 

Other robustified, open-loop approaches, however, attempt to find a smooth contin- 
uous forcing function (e.g., a versine function) that begins and ends with zero slope. 
The basic idea behind such approaches is that a smooth control input without sharp 
transitions is less likely to excite structural modes during maneuvers. On the contrary 
to such a common notion, the results of this paper indicate that properly modulated, 
on-off pulse sequences can achieve a fast maneuvering time with a minimum of struc- 
tural vibrations during and/or after a maneuver, even in the face of plant modeling 
uncertainty. 

The remainder of this paper is organized as follows. Section 2 describes the standard 
time-optimal control problem of flexible spacecraft without modeling uncertainty. A 
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parameter optimization problem is formulated, in which the objective function is the 
maneuvering time. Three cases are explored, as illustrated in Fig. 1, in order to assess the 
actuator placement problem for time-optimal control of flexible spacecraft with multiple 
jets. In Section 3, we investigate the same problem as in Section 2 but considering the 
presence of modeling uncertainty. A time-optimal control problem subject to additional 
robustness (or sensitivity) constraints is formulated, and numerical solutions for three 
cases are then compared with solutions obtained in Section 2. 

2. Time-Optimal Rest-to-Rest Maneuver Control 

Problem Formulation 

Consider a flexible spacecraft described by 

Mx + K x = Gu (1) 

where x is a generalized displacement vector, M a mass matrix, K a stiffness matrix, G 
the control input distribution matrix, and u the control input vector. 

Equation (1) is transformed into the modal equations: 

jil + ^lJ/l = + <t>l2 u 2 + <^13^3 

h + w 2l/2 = ^21 u l + 4>22^2 + <t> 23«3 ^ 

y n + uh/n = 4>n\U\ + d»n2«2 + <An3«3 

where t/< is the i tfl modal coordinate, u>i the i tfl modal frequency, <f>ij the modal input dis- 
tribution coefficient, and n the number of modes considered in control design. Without 
loss of generality, only three control inputs are considered in Eq. (2). 

In this paper, we consider a simple model which is a generic representation of a 
spacecraft with a rigid-body mode and two flexible modes, as shown in Fig. 1. Three 
cases are studied: (i) Case 1 with both “positive” and “negative” jets placed at body 1, 
(ii) Case 2 with a “positive” jet at body 1 and a “negative” jet at body 2, and (iii) Case 

3 with a “positive” jet at body 1 and a “negative” jet at body 3. Case 1 is a typical 

case in which two opposing jets are colocated. In Cases 2 and 3, two opposing jets are 
not colocated. 
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For a class of problems, such as Case 1, the modal equations become 

y. + = fau; i = l , -..,71 ( 3 ) 

where <j> t the i th modal gain and the scalar control input u is bounded as 

- 1 < u < 1 (4) 

The problem is to find the control input u(t) which minimizes the performance index 

J = I dt = t, 

Jo 

subject to Eqs. ( 3 ) and ( 4 ), and given boundary conditions. The time-optimal bang-bang 
solution of this problem with the rest-to-rest maneuvering boundary conditions has, in 
most cases, (2 n - 1) switches, and the solution is symmetric about the mid-maneuver 

time tjj 2 [2]. That is, for a case with (2n - 1) switches, we have the symmetric switching 
pattern given as: 

t i - ij in - t 2n -j; j = 1, • • ,n (5) 

where t } is the j tK switching time and t 2n = t } . 

A bang- bang input with ( 2 n — 1) switches is expressed as 

2n 

u (0 = ££»“•(<-*>) ( 6 ) 

j=o 

where B, is the magnitude of a unit step function u,(t ) at t r This function can be 
characterized by its switch pattern as: 

B={B 0 , B u B 2 , •••, 5 2n } 

T={t 0 , t u t 2 , •••,t 2n } 

where B represents a set of B, with B 0 = B 2n = ±1 and Bj = ±2 for j = 1, • • • , 2 n - 1; 

T re P resents a set of switching times (fj, • • • , t 2n -i) and the initial and final times (t 0 = 0 
and t s = t 2n ). 

The rest-to-rest maneuvering constraint for the rigid-body mode (u>i = 0) can be 
found as 

-!/,(</) = 0 ( 7 ) 

]=0 
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The i tfi structural mode solution for the control input of Eq. (6) for t > tf is 


y,(0 = cosu, <( f ~ *;) 

j=0 

i . 2n 

= j[coso;i(^ — < n )S 5 J cosu; *(^ “ *") 


w; 


j=0 


2n 


+ sinu>i(< - i n ) £ Bj sinu\(fj - *„)] 

3=0 

and it can be shown that the following constraint equation for each mode 

2n 

^2 Bj sin U3i(tj -tn) = 0 

3=0 


( 8 ) 


( 9 ) 


is always satisfied for any bang-bang input which is symmetric about the mid-maneuver 
time t n ■ Consequently, we have the following flexible mode constraints for no-residual 
structural vibration (i.e., yi(t) = 0 for t > tf): 


2 n 


^2 Bj cos u>i(tj - <„) = 0 

3=0 


( 10 ) 


for each flexible mode. 

On the other hand, for Case 2 with the pulse sequences as illustrated in Fig. 2, 
the boundary conditions of the rest-to-rest maneuvering problem result in the following 
constraint: N 

<t>u XI Aj ~ 4>u XI ^3 = 9 ( ^ ) 

3 = 0,2 3 = 1.3 

where <fru and are the modal input distribution coefficients associated with the rigid- 
body mode and the two control inputs iq and u 2 , and Aj and N are defined as in Fig. 2. 

The positioning constraint for the rigid-body mode, with specified y\(t > </), then 
becomes 


jV-l 

2j/l(f/) ~ <?hl XI i 2t 3^3 ~ A j] 

3 = 0,2 

+ 4>\2 X] ~ A»1 = ^ 

3= 1.3 


( 12 ) 


Also, rest-to-rest maneuvering requires that for each flexible mode, yi(t) — 0 for t > t /; 
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i.e., we have 


N - 1 N 

~ <t>n E^ + &2 EE Cij =0 (13a) 

i= o.2 j=l,3 

n-i a/ 

^.1 EE 3 *i “ ^<2 EE 5 »i = 0 (13b) 

J=0,2 j = 1 ,3 

where 


Cjj cos(u>j<j) cos(u; t (<j -|- Aj)) 

5,j = sin(u;^j) — sin(o\(tj + A x )) 

for i = 2, • • • , n. 

Remark: For cases in which the control inputs are one-sided, each control input for 
the time-optimal solution need not be an odd function about the mid-maneuver time. 
Thus, the problem with one-sided control inputs becomes more difficult to solve than 
the standard problem with two-sided control inputs, and many theoretical issues (e.g., 
the uniqueness and structure of time-optimal solutions) need to be resolved. 

We now present a detailed solution of each case. 

Case 1 with a Two-Sided Control Input 

In this case, as illustrated in Fig. 1, a two-sided control input is bounded as: 

-1 < u(t) < +1 

and the equations of motion for this case are 


+ ki(x\ — x 2 ) = u 
TTI 2 X 2 + k\(x 2 — ij) + k 2 (x 2 — x 3 ) = 0 
m 3^3 + k 2 {x^ ~ X 2 ) = 0 

where x x , x 2 and x 3 are the positions of body 1, body 2 and body 3, respectively, and 
the nominal parameters are mi = m 2 = m 3 = *, = k 2 = 1 with appropriate units, and 
time is in units of second. 
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The boundary conditions for a rest-to-rest maneuver are given as 


Zx(0) = *a(0) = X 3 (0) = 0 (14a) 

x^tf) = x 2 {tf) = x 3 {tf) = 1 ( 14b ) 

ii(0) = x 2 (0) = i 3 (0) = 0 (14c) 

Xi(tf) = x 2 (tf) = x 3 (tf) = 0 (14d) 

The modal equations are 

y, = 0.3333u (15a) 

y 2 + W2I/2 = 0.5u (15b) 

y 3 + u>ly 3 = — 0.1667u (15c) 


where 

u\ = [-6 - \fb 2 — c]/2 

u\ = [-b+yJb^4kM/2 

, -^(m, + m 2 )m 3 - ^(m 2 + m 3 )m, 

b — - 

m 1 m 2 m 3 

C = (mi + m 2 + m 3 ) / (m 1 m 2 m 3 ) 
and u>2 = 1 rad/sec and w 3 = \/3 rad/sec for the nominal system. The corresponding 


boundary conditions for the modal coordinates are 

l/i(0) = 2/2(0) = y 3 (0) = 0 (17a) 

y2(tj) = s/3(*/) = 0 ( 17b ) 

J/i(0) = 5/2(0) = y 3 (0) = 0 (17c) 

yi(*/) = ini*}) = &(*/) = 0 ( 17d ) 

The time-optimal control input with 5 switchs is expressed as 

u(£) = u,(t) — 2 u,(t — <i) + 2 u a (t — <2) — 2u,(i — t 3 ) 

+ u,(t — <4) — u,(t — <s) + u s(t ~ te) (1^) 


(16a) 

(16b) 

(16c) 

(16d) 
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with the symmetry conditions 


^6 — 2^3 
t$ = 2t 3 — 1\ 
ti = 2t 3 — 


The time-optimal control problem is then formulated as the following constrained 
minimization problem: 


min J = t / = < 6 = 2t 3 
subject to the following constraints: 


(19) 



( 20 ) 


1 + COs(u? 2^3) + 2 ^(-l)-' cos(u> 2 (< 3 ~ <;)) = 0 

j = 1 
2 

1 + cos(u> 3 t 3 ) + m -1) J cos(u> 3 (* 3 - tj)) = 0 
;=i 

*lt *2, i 3 > 0 

A standard IMSL FORTRAN subroutine was used to obtain a solution as 


( 21 ) 

(22) 

(23) 


u(t) = u s (t) - 2 u 3 (t - 0.944) + 2u s (t - 2.012) 

- 2 u,(t - 3.255) + u 3 (t - 4.499) 

— u s (t — 5.567) + u a (t — 6.511) (24) 

The time responses of x 3 to this time-optimal control input are shown in Fig. 3 for 
four different values of k = k x = k 2 . We notice that the responses are quite sensitive to 
variations in the model parameters. Similar responses can also be observed for arbitrarily 
combined variations of fc, and m,, but keeping the total mass constant (m 1 +m 2 +m 3 = 3). 

For convenience, simulation results only for k = = k 2 variations are presented in this 

paper. 


Case 2 with Two One-Sided Control Inputs 

For Case 2, as illustrated in Fig. 1, two one-sided control inputs are bounded as 

0 < Ul < +1 
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(25a) 



- 1 < U 2 < o ( 25b ) 

For this case, the time-optimal control inputs are assumed as: 

ui = u s (t) — u,(t — Ao) 


+ U,(t - < 2 ) - U,(t -ti~ A 2 ) 
u 2 = — u,(t — ti) + u s (t — t\ — Ai) 

(26a) 

— u,(t — t 3 ) + u,(£ — t 3 — A 3 ) 
The modal equations for this case are 

(26b) 

y x = 0.3333UJ + 0.3333u 2 

(27a) 

j ?2 + W 2 J /2 = 0.5ui 

(27b) 

y 3 + w|y 3 = 0.1667u, - 0.3333 u 2 
and the rest-to-rest maneuver constraints are 

(27c) 

Ao — Ai + A 2 — A 3 = 0 

(28a) 

6 + D-l)iA?- 2 ( ) A J ! = 0 

7 = 0 

(28b) 

J — ^ 

52 [cos(w 2 fj) - cos(w 2 (ij + Aj))] = 0 
3- 1.3 

(28c) 

52 (sin(u; 2 tj) — sin(u; 2 (tj 4- Aj))] = 0 

j = l,3 

52 [cos(u> 3 tj) - cos(u> 3 (t_, + Aj))] 
j=1.3 

(2Sd) 

+ 2 51 [cos(u } 3 tj) - cos(a) 3 (tj + Aj))] = 0 
i= 0,2 

52 [sin(w 3 <j) — sin(w 3 (<j + Aj))] 
j= 1.3 

(28e) 

+ 2 5Z [ s > n ( w 3 <j) — sin(aj 3 (<j + Aj))] = 0 

j= 0,2 

(28f) 

The time-optimal solution can then be obtained by solving the constrained minimization 
problem: 

min J = tj = t 3 + A 3 

subject to the constraint given by Eqs. (28). 

(29) 
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A solution of this problem can be found as 


t 0 = 0.0000, A 0 = 0.8459 
= 1.1594, A, = 1.0487 

t 2 = 2.1455, A 2 = 1.2516 

<3 = 4.3010, A 3 = 1.0487 
t f = 5.3497 


(30) 


The time responses of x 3 to the time-optimal control inputs are shown in Fig. 4 for 
four different values of k — k\ — k 2 . Similar to Case 1, the responses are sensitive to 
variations in the model parameters. An interesting feature of this case is that the pulse 
sequences are of a “bang-off-bang” type, resulting in the control on-time of 4.195 sec, 
which is different from the maneuver time of 5.35 sec. 


Case 3 with Two One-Sided Control Inputs 

For Case 3, as illustrated in Fig. 1, two one-sided control inputs are bounded as 


0 < ui < +1 
— 1 < u 3 < 0 

Similar to Case 2, the time-optimal control inputs are assumed as: 


(31a) 

(31b) 


Ui = u.(<) - u,(t - A 0 ) 

+ «•(<- *2) -u.(t-t 2 - A 2 ) ( 32 a) 

«3 = - u,(t ~ <i) + u 4 (< - t, - AO 

— u a (t — f 3 ) + u s (t — t 3 — A 3 ) (32b) 

The modal equations for this case are 

yi = 0.3333(ui 4- u 3 ) (33a) 

£2 + u>ly 2 = 0.5(«i - u 3 ) (33b) 

V3 + w 3 j/ 3 = 0.1667(u! + u 3 ) (33c) 

and the rest-to-rest maneuver constraints are 

Aq — A) 4 A 2 — A 3 = 0 (34a) 
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(34b) 


3 = 0 


3 


^( — 1) J [sin(u> 2 <j) - s\n(u3 2 (tj -1- A,))] = 0 

3=0 

(34c) 

52(-iy[cos(u; 2 <j) - cos{u 2 {tj + Aj))] = 0 

3=0 

(34d) 

E( — 1) ; [sin(uj 3 <j ) - sin(u> 3 (<_, + Ay))] = 0 

3=0 

(34e) 

^(-l) J [cos(u> 3 fy) - cos(u3 3 (tj + Ay))] = 0 

( 34 f ) 


3=0 


tj, hi ^3) ^4) ^5 >0; ^0 — 0 

A s > 0 


The time-optimal solution can then be obtained by solving the constrained mini- 
mization problem: 

min J = tj = < 3 + A 3 (35) 

subject to the constraints given by Eqs. (34). 

A solution of this problem is 

t 0 = 0.0000, A 0 = 0.9510 

t\ = 1.3631, A! =0.1658 

t 2 = 2.8329, A 2 = 0.1658 (36) 

< 3 = 3.4109, A 3 = 0.9510 
t s = 4.3619 

The time responses of x 3 to the time-optimal control inputs are shown in Fig. 5 for 
four different values oi k = k x = k 2 . Again, the responses are quite sensitive to variations 
in the model parameters. Similar to Case 2, an interesting feature of this case is that 
the pulse sequences are of a “bang-off- bang” type, resulting in the control on-time of 
2.234 sec, which is different from the maneuver time of 4.362 sec. 

Compared to Case 1 and Case 2, this case has the fastest maneuver time as well 
as the smallest control on-time. Therefore, the actuator configuration of Case 3 can be 
considered to be “optimal” in the sense of minimizing both the maneuver time and the 
jet on-time. 
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3. Robust Time-Optimal Control 

As shown in the preceding section, a standard, time-optimal control approach re- 
quires an accurate mathematical model and thus the resulting solution is often sensitive 
to plant modeling uncertainty. 

In this section, expanding on the approach introduced in Section 2, a parameter 
optimization problem is formulated with additional constraints for robustness with re- 
spect to the structural frequency uncertainty. The resulting robustified or desensitized , 
time-optimal soultion is a multi-switch bang-bang control, and is thus implementable 
for spacecraft equipped with on-off reaction jets [3]. 


Problem Formulation 

By taking the derivative of Eq. (8) with respect to w,, we get 


dyi(t) <t>> 


duii uf C ° S ~2 ^ ^ tj ~ ~2 ^ Bj sin 

for each flexible mode. Letting rfy, (<)/<&;, = 0 for all t > t h we have 


2n 


</> 


t f^ 


(37) 


2 n 


E(*> - %) B 1 cos Ufa - l -L) = 0; i = 2, • • • , n 


tf 


>=0 


(38) 


which are called the first-order robustness constraints for the case with a two-sided 
control input. 

Similarly, the robustness constraints for a case with two one-sided control inputs 


can 


be found as: 


Af-l N 

- <t>i\ E c 'j + ta E c 'j = 0 

J=0,2 >=1.3 

N - 1 N 

<t>n s 'j - ^>2 E = o 

]= 0,2 >= 1,3 


(39a) 

(39b) 


where 


Cij = *> cos(u >itj) - (tj + A,) cos (u>i(tj + Aj)) 
s i} = tjsiniuitj) - (tj + A i )sin(a; l (< i + A,)) 
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which are called the first-order robustness constraints for the case with two one-sided 
control inputs. 

Case 1 with a Two-Sided Control Input 

For Case 1, the time-optimal control is a five-switch bang-bang function, but the 
resulting responses were shown to be very sensitive to variations in model parameters. 
A robustified, time-optimal solution of the same problem is now computed as follows. 
The robust time-optimal control input is assumed as: 

u(t) = u,(<) + 2 — ^i)l d" u,(< — <io) (40) 

j=i 

with the symmetry conditions 

<6 = 2<5 — 1 4 
tj = 2< 5 "" ^3 

<8 — 2< 5 — h 

<9 = 2<5 — <1 
<io = 2< 5 

The constrained optimization problem with the first-order robustness constraint can 
be formulated as: 

min J = tf = t io (42) 

subject to 

6 + E(— l) J+, [2(<io-^) 2 ]-4 = 0 

;=l 

1 + cos(u; 2 <s) + 2£[(-l) J cos (u 2 (tj ~ h))] = 0 

j=i 

4 

1 + cos(u> 3 < 5 ) + 2^[(-l) J cos(u) 3 (< j — < 5 ))] = 0 
}= 1 

< 5 sin(u>2<5) + 2 — <5) sin(ui2(<j — <5)) = 0 

j = 1 

< 5 sin(u> 3 < 5 ) + 2^(-l) j (C - < 5 ) sin(oj 3 (<j - < 5 )) = 0 

J = 1 

< o = 0 
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A robust time-optimal solution with 9 switches can be found as: 


<1 = 0.560, * 2 = 1.460 
< 3 = 2.690, < 4 = 3.804 

*5 = 5.091, t 6 = 6.377 (43) 

*7 = 7.491, *8 = 8.722 
*9 = 9.622, * 10 = 10.18 

The time responses of x 3 to this robust time-optimal control input are shown in Fig. 6 
for four different values of k = k x = k 2 . We notice that the resulting responses are less 
sensitive to parameter variations, as compared to the responses to the ideal, time-optimal 
control input, as shown in Fig. 3. The second flexible mode is significantly excited during 
manuevers, however. Performance robustness has been increased at the expense of the 
increased maneuvering time of 10.18 sec, as compared to the ideal minimum-time of 6.511 
sec. It is, however, emphasized that simply prolonging the maneuver time does not help 
to reduce residual structural vibrations caused by modeling uncertainty; a proper pulse 
sequence is necessary. 

Case 2 with Two One-Sided Control Inputs 

For Case 2, we can represent the control inputs as follows: 

ui = «,(*) - u,(* - A 0 ) -I- u s (t - * 2 ) 

- u,(* - < 2 - A 2 ) + u,(< - * 4 ) - U,(f - t 4 - A 4 ) 

U 2 = - u s (t - *j) + u s (t - *, - A,) - u s (< - *3) 

+ U»(* - *3 - A 3 ) - u,(t - * 5 ) + U,(t - *5 - As) 

where we have 11 unknowns to be determined, and t } and A, are defined as in Fig. 2. 

The robust time-optimal solution can then be obtained by solving the constrained 
parameter optimization problem 


min J = tf = * 5 A s 


( 44 ). 


subject to 


Aq — Ai + A 2 — A 3 + A 4 — A5 = 0 
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6 + 2<jA ; ] — 0 

>=o 

y [cOs(u>2<j) — COs(u;2(<j + Aj))] =0 
>=1.3 

52 [sin(u> 2 f>) - sin(u> 2 (<> + Aj))] = 0 

>=1.3 

5^ [oos(u>3<j) - COs(w 3 (<j + A ; ))] 

>=1.3 

+ 2 51 [cO s ( w 3^) - COs(u> 3 (<j + Aj))] = 0 

j=0,2,4 

52 [sin(a!3<j) -sin(w 3 (<j + Aj))] 

>=1,3 

+ 2 5Z (sin(u> 3 *>) -sin(u> 3 (<j + Aj))] = 0 

>= 0 , 2,4 

52 [*> COs(u > 2 <>) — (tj + Aj) COs(w 2 (<j + Aj))] = 0 

>=1,3 

52 [<> sin(u; 2 <>) - (tj + Aj) sin(u; 2 (<, + Aj))] = 0 

>=1,3 

52 [tj cos(u;3 tj) - (t } + Aj) cos(u> 3 (^ + Aj))] 

>=1,3 

+ 2 52 [t jCOs(u>3*j) - (<; + A ; )cos(u;3(<j + Aj))] = 0 

>=0,2,4 

52 [tj sin(u>3<j) - {tj + Aj) sin(u? 3 (fj + Aj))] 

>=1,3 

+ 2 52 (*>sin(w3<>) - (<> + Aj) sin(w 3 (<j + A ; ))] = 0 

>=0,2,4 

Aj > 0; j = 0,1, 2, 3, 4,5 

< 1 , t 2 , <3, U, ts > 0 ; <0 = o 

A solution to this problem is: 

<o = 0.0000, A 0 = 0.3855 
*, = 0.4440, Ai = 0.5783 
< 2 = 1.6118, A 2 = 1-4181 

<3 = 3.2990, A 3 = 1-2131 (45) 

U = 4.8208, A 4 = 0.5661 
< 5 = 6.7887, A 5 = 0.5783 
tj = 7.3671 

The time responses of 13 to the robust time-optimal control inputs are shown in 
Fig. 7 for four different values of k = k x = k 2 . Similar to Case 1 of the preceding section, 
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the robustness has been increased at the expense of the increased maneuvering time of 
7.367 sec, as compared to the ideal minimum-time of 5.35 sec. The jet on-time is 4.195 
sec. The second flexible mode is less excited, as compared to Case 1. 

Case 3 with Two One-Sided Control Inputs 

Assuming that each control input has two pulses as in Case 2, we can represent the 
control inputs as: 


«i = MO - M< “ A 0 ) + M* - t 3 ) - u a (t -t 2 - A 2 ) 

+ M* — t 4 ) — u a (t - t 4 — A 4 ) 
u 3 = - u a (t - tj) + M< - 0 ~ A,) - u,(t - t 3 ) 

-I- u a {t - <3 - A 3 ) - U a (t - t s ) + U a (t -t 5 - As) 

where we have 11 unknowns to be determined, and t, and A; are defined as in Fig. 2. 

The robust time-optimal solution can then be obtained by solving the constrained 
parameter optimization problem 

min J = t s + A s (47) 

subject to 


Ao — Ai + A2 — A 3 + A 4 — A5 = 0 
6 + D- 1 ) > [A j 2 -2CA;] = 0 

j=o 

5 

- sin(u> 2 (0 + Aj))] = 0 

j=o 

5 

^ ( — 1 ) J [cos (u> 2 tj ) - cos(wj(tj + Aj))] = 0 

j=o 

5 

^[(-1)-' sin (u 3 tj) - sin(u; 3 (C + Aj))] = 0 

j=o 

5 

cos(w 3 C) - cos (w 3 (C -1- Aj))] = 0 

i=o 

5 

^(-ly^coswatj - (t } + A : )cosu 2 (t : -f A,)] = 0 
]=0 
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y^(— iyfo sinu >2 tj — {tj + A ; ) s\nui 2 {t } + Aj)] — 0 

1=0 

£(-1)-% cosu> 3 tj - {tj + Aj) cos w 3 (fj + Aj)] = 0 

i=o 

£(-l) J [fj sin u 3 tj - {tj + Aj) sin u 3 (tj + Aj)] = 0 

i=o 

Aj > 0; j= 0,1, 2, 3, 4, 5 

t\y ^2) ^3) ^4i ^5^0, <0 0 

A solution to this problem is 

< 0 = 0.000, A 0 = 0.2189 
h = 3.774, A x = 0.2609 
t 2 = 2.030, A 2 = 0.3594 

t 3 = 5.778, A 3 = 0.3594 (48) 

t A = 4.152, A< = 0.2609 
* 5 = 7.968, A 5 = 0.2189 
t j — 187 

The time responses of x 3 to the robust time-optimal control inputs are shown in Fig. S 
for four different values of it = k, = k 2 . The robustness has been increased at the expense 
of the increased maneuvering time of 8.187 sec, as compared to the ideal minimum- 
time of 4.362 sec. However, the control on-time is only 1.678 seconds! Because of the 
properly coordinated pulse sequences, the flexible modes are not significantly excited 
during manuevers and the residual responses after the maneuvers are well desensitized. 


4. Summary 

In this paper, we have demonstrated that the proposed robustification or desensitiza- 
tion approach does generate robust time-optimal open-loop control inputs for uncertain 
dynamical systems. Furthermore, on the contrary to a common notion, the results of 
this paper indicate that properly coordinated, on-off pulse sequences can achieve a fast 
maneuvering time with a minimum of structural vibrations during and/or after a ma 
neuver, even in the face of plant modeling uncertainty. The time-optimal responses 
have been desensitized at the expense of the increased maneuvering time. It is again 
emphasized that simply prolonging the maneuver time does not help to reduce residual 
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structural vibrations caused by modeling uncertainty; a proper coordination of pulse 
sequences is necessary, as demonstrated in this paper. 

The results of this paper are summarized in Table 1. As can be noticed in this 
table, it is natural to select the actuator configuration of Case 3, since this case provides 
the “best” overall performance in the sense of minimizing the maneuvering time, fuel 
consumption (jet on-time), and structural mode excitation. For Case 1, the maneuvering 
time and the jet on-time are the same, which is clearly undesirable from the viewpoint of 
fuel consumption. To avoid such undesirable continuous jet firings during a maneuver, 
a robust fuel- and time-optimal control problem is formulated in [4]. 

5. Conclusions 

A time-optimal open-loop control problem of flexible spacecraft in the face of mod- 
eling uncertainty has been investigated. The primary study objective was to explore 
the feasibility of computing open-loop, on-off pulse control logic for uncertain flexible 
spacecraft. The results indicate that the proposed approach significantly reduces the 
residual structural vibrations caused by modeling uncertainty. The results also indicate 
the importance of a proper jet placement for practical trade-ofTs among the maneuvering 
time, fuel consumption, and performance robustness. 
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